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Abstract 

The Vlasov equation for the coUisionless evolution of the single-particle probability distri- 
bution function (PDF) is a well-known Lie-Poisson Hamiltonian system. Remarkably, the 
operation of taking the moments of the Vlasov PDF preserves the Lie-Poisson structure. The 
individual particle motions correspond to singular solutions of the Vlasov equation. The pa- 
per focuses on singular solutions of the problem of geodesic motion of the Vlasov moments. 
These singular solutions recover geodesic motion of the individual particles. 
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1 Introduction 

The Vlasov equation. The evolution of N identical particles in phase space with coordinates 
{qi,Pi) i = 1,2, . . . ,N , may be described by an evolution equation for their joint probability 
distribution function. Integrating over all but one of the particle phase-space coordinates yields 
an evolution equation for the single-particle probability distribution function (PDF). This is the 
Vlasov equation, which may be expressed as an advection equation for the phase-space density / 
along the Hamiltonian vector field X/^ corresponding to single-particle motion with Hamiltonian 
H{q,p): 

^={f,H}= -diV(,,,)(/X^) =-C^,f with XHiq,p) = (^,-^) (1-1) 

The solutions of the Vlasov equation reflect its heritage in particle dynamics, which may be 
reclaimed by writing its many-particle PDF as a product of delta functions in phase space. Any 
number of these delta functions may be integrated out until all that remains is the dynamics of 
a single particle in the collective field of the others. 

In the mean-field approximation of plasma dynamics, this collective field generates the total 
electromagnetic properties and the self-consistent equations obeyed by the single particle PDF 
are the Vlasov-Maxwell equations. In the electrostatic approximation, these reduce to the 
Vlasov-Poisson (VP) equations, which govern the statistical distributions of particle systems 
ranging from integrated circuits (MOSFETS, metal-oxide semiconductor field-effect transistors), 
to charged-particle beams, to the distribution of stars in a galaxy. 

A class of singular solutions of the VP equations called the "cold plasma" solutions have a 
particularly beautiful experimental realization in the Malmb erg-Penning trap. In this experi- 
ment, the time average of the vertical motion closely parallels the Euler fluid equations. In fact, 
the cold plasma singular Vlasov-Poisson solution turns out to obey the equations of point-vortex 
dynamics in an incompressible ideal flow. This coincidence allows the discrete arrays of "vortex 
crystals" envisioned by J. J. Thomson for fluid vortices to be realized experimentally as solutions 
of the Vlasov-Poisson equations. For a survey of these experimental cold-plasma results see (ll] . 

Vlasov moments. The Euler fluid equations arise by imposing a closure relation on the flrst 
three momentum moments, or p— moments of the Vlasov PDF f{p, q, t). The zero-th p— moment 
is the spatial density of particles. The first p— moment is the mean fiuid momentum. Introducing 
an expression for the fluid pressure in terms of the density and momentum closes the system of 
moment equations, which otherwise would possess a countably inflnite number of dependent 
variables. 

The operation of taking p— moments preserves the geometric nature of Vlasov dynamics. 
In particular, this operation is a Poisson map. That is, it takes the Lie-Poisson structure for 
Vlasov dynamics into another Lie-Poisson system. However, strictly speaking, the solutions for 
the p— moments cannot yet be claimed to undergo coadjoint motion, as in the case of the Vlasov 
PDF solutions, because the group action underlying the Lie-Poisson structure found for the 
p— moments is not yet understood. 

A closure after the first p— moment results in Euler's useful and beautiful theory of ideal 
fluids which is also Lie-Poisson. As its primary geometric characteristic, Euler's fiuid theory 
represents fluid flow as Hamiltonian geodesic motion on the space of smooth invertible maps 
acting on the flow domain and possessing smooth inverses. The (left) action of these smooth 
maps (called diffeomorphisms) on the fiuid reference configuration moves the fiuid particles 
around in their container. And their smooth inverses recall the initial reference configuration 
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(or label) for the fluid particle currently occupying any given position in space. Thus, the motion 
of all the fluid particles in a container is represented as a time-dependent curve in the infinite- 
dimensional group of diffeomorphisms. Moreover, this curve describing the sequential actions 
of the diffeomorphisms on the fluid domain is a special optimal curve that distills the fluid 
motion into a single statement. Namely, "A fluid moves to get out of its own way as efficiently 
as possible". Put more mathematically, fluid flow occurs along a curve in the diffeomorphism 
group which is a geodesic with respect to the metric on its tangent space supplied by its kinetic 
energy. 

Given the beauty and utility of the solution behavior for Euler's equation for the first 
p— moment, one is intrigued to know more about the dynamics of the other moments of Vlasov's 
equation. Of course, the dynamics of the the p— moments of the Vlasov-Poisson equation is one 
of the mainstream subjects of plasma physics and space physics. 

Summary. This paper formulates the dynamics of Vlasov p— moments governed by quadratic 
Hamiltonians. This dynamics is a certain type of geodesic motion on the symplectomorphisms, 
rather than on the diffeomorphisms for fluids. The symplectomorphisms are smooth invertible 
maps acting on the phase space and possessing smooth inverses. The theory of moment dynam- 
ics for the Vlasov equation turns out to be equivalent to the theory of shallow water equations, 
and a particular example is the one-dimensional system of Benney long wave equations, which 
is integrable [U [12] . 

Here we shall consider the singular solutions of the geodesic dynamics of the Vlasov p— moments. 
Remarkably, these equations turn out to be related to other integrable systems governing shal- 
low water wave theory. For example, when the Vlasov p— moment equations for geodesic motion 
on the symplectomorphisms are closed at the level of the first p— moment, their singular solu- 
tions are found to recover the peaked soliton of the integrable Camassa-Holm (CH) equation for 
shallow water waves [5]. These singular Vlasov moment solutions also correspond to individual 
particle motion. 

Thus, geodesic symplectic dynamics of the Vlasov p— moments is found to possess singular 
solutions whose closure at the fluid level for the CH equation recovers the peakon solutions of 
shallow water theory. Being solitons, the CH peakons superpose and undergo elastic collisions in 
fully nonlinear interactions. The singular solutions for Vlasov p— moments presented here also 
superpose and interact nonlinearly as coherent structures. 

The plan of the paper follows: 

Section [2] reviews the Vlasov p— moment equations and recounts their Lie-Poisson Hamiltonian 
structure using the Kupershmidt-Manin Lie-Poisson bracket. Variational formulations of 
the p— moment dynamics are also provided. 

Section [3] shows how the Lagrangian framework for fluid dynamics is recovered from the Vlasov 
p— moments and establishes connections with some equations for shallow water waves. In 
this case, the p— moment equations are shown to possess singular solutions. 

Section [4] establishes the connections between the integrable Benney equations and the physics 
of charged-particle accelerator beams. To our knowledge, these connections are noted here 
for the first time. We also point out how the experimental realization of solitary waves in 
coasting particle beams has its roots in the integrability of the Benney system. 

Section [5] formulates the problem of geodesic motion on the symplectomorphisms in terms of 
the Vlasov p— moments and identifies the singular solutions of this problem. This geodesic 
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motion is related to a geodesic form of the Vlasov equation. Thus the singular solutions are 
found to originate in the single particle dynamics on phase space. In a special case, the 
truncation of geodesic symplectic motion to geodesic diffeomorphic motion for the first 
p— moment recovers the singular solutions of the Camassa-Holm equation, and thereby 
correspond to single particle dynamics. 

The geodesic form of the Vlasov equation was introduced in [15j, where it was also shown 
how to extend the treatment to higher dimensions. 



2 Review of Vlasov moment dynamics 



The Vlasov equation may be expressed as 



dl 
dt 



f 



ad 



■5h/5f 



f 



(2.1) 



5hi _df d 6h df d 6h 
6f J dp dq 6f dq dp 6f 

Here the canonical Poisson bracket [ • , • ] is defined for smooth functions on phase space with 
coordinates {q,p) and f{q,p,t) is the evolving Vlasov single-particle distribution function. The 
variational derivative Sh/6f is the single particle Hamiltonian and the adsh/Sf f explained as 
follows. 

A functional g[f] of the Vlasov distribution / evolves according to 



dg^ 
dt 



/ 



6g 6h 
If' If 



dqdp 



6h 
/, 



dqdp 

~6g 6h 
IJf'Sf 



{9,h} 



In this calculation boundary terms were neglected upon integrating by parts in the third step 
and the notation (( • , • )) is introduced for the pairing in phase space. The quantity {g , h} 
defined in terms of this pairing is the Lie-Poisson Vlasov (LPV) bracket [3U] . This Hamiltonian 
evolution equation may also be expressed as 



'^9 f > T 



f , &dsh/5f 



6f// 



ad* sh/5f f , 



Sg_\s 
6f// 



which defines the Lie-algebraic operations ad and ad* in this case in terms of the pairing on 
phase space (( ■ , • )) : 5* x s i— 



The notation adj^^^j / in (j2.ip expresses coadjoint action of 
Sh/5f G 5 on / G s*, where s is the Lie algebra of single particle Hamiltonian vector fields and 
s* is its dual under pairing in phase space. This is the sense in which the Vlasov equation 
represents coadjoint motion on the symplectomorphisms. This Lie-Poisson structure has also 
been extended to include Yang-Mills theories in [13] and |14j . 

In higher dimensions, particularly n = 3, we may take the direct sum of the Vlasov Lie- 
Poisson bracket, together with with the Poisson bracket for an electromagnetic field (in the 
Coulomb gauge) where the electric field E and magnetic vector potential A are canonically 
conjugate. For discussions of the Vlasov-Maxwell equations from a geometric viewpoint in the 
same spirit as the present approach, see [7], [25], [26] and [30] . 

2.1 Dynamics of Vlasov Moments 

The phase space moments of the Vlasov distribution function are defined by 



9 mm 



= 11 f{q,p)q'^p"'dqdp. 
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The moments gfhm ^-re often used in treating the colhsionless dynamics of plasmas and 
particle beams [10]. This is usually done by considering low-order truncations of the potentially 
infinite sum over phase space moments, 

oo oo 

9 ^ ^ ^ fhm9 rh,m ) ^ ^ ^ ^ nn9 n,n i 

rh,m=0 n,n=0 

with constants flmm ^nm with m.,m,n,n = 0,1, ... . If /i is the Hamiltonian, the sum over 
moments g evolves under the Vlasov dynamics according to the Poisson bracket relation 

dg 

"^ = {5)^} = (^mmbnn{^^~^'IT')9fh+n-l,rn+n-l- 
rh,m,n,n=0 

The symplectic invariants associated with Hamiltonian flows of the moments were dis- 
covered and classified in [18j . Finite dimensional approximations of the whole q,p— moment 
hierarchy were discussed in For discussions of the Lie-algebraic approach to the control 

and steering of charged particle beams, see [TO]. 



2.2 Dynamics of Vlasov p— Moments 

In contrast to the (7, p— moments, the momentum moments, or "p— moments," of the Vlasov 
function are defined as 



A^{q,t)= J p"" fiq,p,t)dp, m = 0,l. 



That is, the p— moments are dependent integrals over p of the product of powers p"*, m = 
0, 1, . . . , times the Vlasov solution f{q,p, t). We shall consider functionals of these p— moments 
defined by, 

00 „^ 00 „ 00 

9 = ^ am{q) f dqdp = ^ I am{q) Am{q) dq =: ^ (Ajn , oi 

m=0 m=0 m=0 

00 00 n 00 

h = Mq)p''fdqdp = Y^ / Pn{q)An{q)dq=:Y,[An,f3n), 

n=0 •^•^ n=0 n=0 

where ( ■ , ■ ) is the pairing on position space. 

The functions Om and /3„ with m,n = 0,1,... are assumed to be suitably smooth and 
integrable against the Vlasov p— moments. To assure these properties, one may relate the 
p— moments to the previous sums of Vlasov moments by choosing 

00 00 
am{q) = ^affimq'^ and /3n(g) = ^ 6nng" • 

m=0 n=0 

For these choices of oim{q) and Pniq), the sums of p— moments will recover the full set of Vlasov 
moments. Thus, as long as the moments of the distribution f{q,p) continue to exist 
under the Vlasov evolution, (q) and /?„ (q) are smooth 

functions whose Taylor series expands the p— moments in the g,p— moments. These functions 
are dual to the p— moments Am{q) with m = 0, 1, . . . under the pairing (• , •) in the spatial 
variable q. In what follows we will assume homogeneous boundary conditions. This means, for 
example, that we will ignore boundary terms arising from integrations by parts. 
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The Poisson bracket among the p— moments is obtained from the LPV bracket may be 
expressed as 



m,n=0 



{g,h}{{A}) = - 2^ J jj-^nj^—Am+n-i + [m + n)Am+n-i 



dq 5Ar, 



dq 



This is the Kupershmidt-Manin Lie- Poisson (KMLP) bracket [22], which is defined for functions 
on the dual of the Lie algebra with bracket 

This Lie algebra bracket inherits the Jacobi identity from its definition in terms of the canonical 
Hamiltonian vector fields. Thus, we have shown the 

Theorem 2.1 (Gibbons |12] ) The operation of taking p— moments of Vlasov solutions is 
a Poisson map. It takes the LPV bracket describing the evolution of f{q,p) into the KMLP 
bracket, describing the evolution of the p— moments An{x). 



A result related to this, for the Benney hierarchy [T], was also noted by Lebedev and Manin 
The evolution of a particular p— moment Am{q,t) is obtained from the KMLP bracket by 

dAm f . , . f d d \ 5h ^ 6h 

^ = { , M = - X. (n^^rn+n-l + mA^^^.^-) ^ = Am , A^} 

These moment equations can also be derived from variational principles, as shown in [T5 
within the Hamilton- Poincare framework [8]. 



3 Moments and cotangent lifts of diffeomorphisms 

As explained in the introduction, a first order closure of the moment hierarchy leads to the 
equations of ideal fluid dynamics. Such equations represent coadjoint motion with respect to 
the Lie group of smooth invertible maps (diffeomorphisms). This coadjoint evolution may be 
interpreted in terms of Lagrangian variables, which are invariant under the action of diffeomor- 
phisms. In this section we investigate how the entire moment hierarchy may be expressed in 
terms of the fluid quantities evolving under the diffeomorphisms and express the conservation 
laws in this case. 



3.1 Lagrangian variables and cotangent lifts 

In order to look for Lagrangian variables, we consider the geometric interpretation of the mo- 
ments, regarded as fiber integrals on the cotangent bundle T*Q of some configuration manifold 
Q. A p— moment is defined as a fiber integral; that is, an integral on the single fiber T*Q with 
base point q & Q kept fixed 

An{q)= [ P''f{q,p)dp 

A similar approach is followed for gyrokinetics in [28j. Now, the problem is that in general the 
integrand does not stay on one same fiber under the action of canonical transformations, i.e. 
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symplectomorphisms are not fiber-preserving in the general case. However, one may avoid this 
problem by restricting to a subgroup of these canonical transformations whose action is fiber 
preserving. The transformations in this subgroup (indicated with T*Diff(Q)) are called point 
transformations or cotangent lifts of diffeomorphisms and they arise from diffeomorphisms on 
points in configuration space [24j, such that 

Qt = qt{qo) 

The fiber preserving nature of cotangent lifts is expressed by the preservation of the canonical 
one-form: 

pt dqt = pt {Qo , Pa ) dqt (^o ) = Po d% 

This fact also reflects in the particular form assumed by the generating functions of cotangent 
lifts, which are linear in the momentum coordinate, i.e. 

d 

H{q,p) = P{q) -^-i pdq = p (3{q) . 

Restricting to cotangent lifts represents a limitation in comparison with considering the whole 
symplectic group. However, this is a natural way of recovering the fluid equations, starting from 
the full moment dynamics. 



3.2 Characteristic equations for the moments 

Once one restricts to cotangent lifts, Lagrangian moment variables may be defined and con- 
servation laws may be found, as in the context of fluid dynamics. The key idea is to use the 
preservation of the canonical one-form for constructing invariant quantities. Indeed one may 
take n times the tensor product of the canonical one-form with itself and write: 

P?{dqtr=pS{dq^r 
One then considers the preservation of the Vlasov density 

ft{quPt) dqt A dpt = /o(go,Po) dq^ A dpo 

and writes 

Pt ft{qt,Pt) {dqtT ® dqt A dpt = Po fo{qo,Po) {dqoT dqo A dpo 
Integration over the canonical particle momenta yields the following characteristic equations 



d_ 
di 



1 dH 
4*^ (qt) {dqtT ® dqt\ = along qt = ^ = P{q) (3.1) 



which recover the well known conservations for fluid density and momentum (n = 0, 1) and can 
be equivalently written in terms of the Lie-Poisson equations arising from the KMLP bracket. 
Indeed, if the vector fleld /3 is identified with the Lie algebra variable (5 = (5i = 5h/5Ai and 
h{Ai) is the moment Hamiltonian, then the KMLP form of the moment equations is 

An + ad*p^An = 0. 

In this case, the KM ad^^ operation coincides with the Lie derivation Hence, one may also 
write it equivalently as 

An + £pAn = 0. 
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This equation is reminiscent of the so called h- equation introduced in ^20j, for which the vector 
field (3 is nonlocal and may be taken as f3{q) = G * An (for any n), where G is the Green's 
function of the Helmholtz operator and the star denotes convolution. When the vector field /? 
is sufficiently smooth, this equation is known to possess singular solutions of the form 

K 

An{q,t) =^Pn,i{t)5{q-Q,{t)) 
i=l 

where the i—th position Qi and weight j of the singular solution for the n— th moment satisfy 
the following equations 

^9 <i=Q. 

Interestingly, for n = 1 (with P{q) = G*Ai), these equations recover the pulson solutions of the 
Camassa-Holm equation, which play an important role in the following discussion. Moreover the 
particular case n = 1 represents the single particle solution of the Vlasov equation. However, 
when n 7^ 1 the interpretation of these solutions as single-particle motion requires the particular 
choice Pn,i = (-Pi)"- For this choice, the n— th weight is identified with the n— th power of the 
particle momentum. 

The equations presented in this section provide a geometric interpretation of the moments in 
terms of the covariant tensor power densities in (j3.ip . When one restricts to cotangent lifts, one 
can specify an action of the group of diffeomorphisms on the moments. However, this action is 
not yet understood at the group level in the general case. The general geometric nature of the Lie 
algebra would be expressed in terms of contravariant tensors dual to the covariant tensor power 
densities. For the cotangent lifts, these contravariant tensors reduce to contravariant vector 
fields. In this case, one is able to characterize their action on the moments as Lie derivatives. 



KMLP bracket and cotangent lifts. We have seen that restricting to cotangent lifts leads 
to a Lagrangian fluid-like formulation of the dynamics of the resulting p-moments. In this case, 
the moment equations are given by the KMLP bracket when the Hamiltonian depends only on 
the first moment (/3i = 5h/5Ai) 



An 





] 


1m„ 


Ml J 



If one now restricts the bracket to functionals of only the first moments, one may check that the 
KMLP bracket yields the well known Lie-Poisson bracket on the group of diffeomorphisms 



6h d 5g d 5h 



6Ai dqSAi 5Ai dqSAi 

This is a very natural step since diffeomorphisms and their cotangent lifts are isomorphic. In 
fact, this is the bracket used for ideal incompressible fluids as well as for the construction of the 
EPDiff equation, which will be discussed later as an application of moment dynamics. 
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4 Applications of the KMLP bracket and quadratic terms 
4.1 The Benney equations 

The KMLP bracket (j2.2|) was first derived in the context of Benney long waves, whose Hamil- 
tonian is 

H = l j{A2{q)+gAl{q))dq. 
The Hamiltonian form dtA^ = {An, H} with the KMLP bracket leads to the moment equations 

dAn , dAn+l , . dAo ^ 

at oq oq 

derived by Benney pTj as a description of long waves on a shallow perfect fluid, with a free surface 
aX y = h{q,t). In his interpretation, the An were vertical moments of the horizontal component 
of the velocity p[q, y, t): 

An= f'{q,y,t)dy. 



The corresponding system of evolution equations for p{q, y, t) and h{q, t) is related by hodograph 
transformation, y = f{q,p',t) dp', to the Vlasov equation 

^+p^-g^^-0 (4 1) 

dt dq dq dp 

The most important fact about the Benney hierarchy is that it is completely integrable. This 
fact emerges from the following observation. Upon defining a function \{q,p,t) by the principal 
value integral, 

X{q,p,t)=p + P —dp, 

J-oo P-P 

it is straightforward to verify ^23j that 

d\ dX_ dAp dX _ ^ 
dt ^ dq ^ dq dp ' 

so that / and A are advected along the same characteristics. 



Applications to coasting accelerator beams. Interestingly, the Vlasov equation ()4.ip re- 
sulting from the hodograph transformation of the Benney equation is exactly the same as the 
equation that regulates coasting proton beams in particle accelerators (see for example p2] 
where a bunching term is also included). 

Now, this is a very remarkable fact: the integrability of the Vlasov-Benney equation implies 
soliton solutions, indications for which seem to have been found experimentally at CERN |21j . 
BNL [2], LANL [9] and FermiLab [27\. (In the last case solitons are shown to appear even when 
a bunching force is present.) These solutions have attracted the attention of the accelerator 
community and considerable analytical work has been carried out over the last decade (see for 
example |31]). Nevertheless to our knowledge, the existence of solitons in coasting proton beams 
has never been related to the integrability of the governing Vlasov equation via its connection to 
the Benney hierarchy,, although this would explain very naturally why robust coherent structures 
are seen in these experiments as fully nonlinear excitations. We plan to pursue this direction in 
future research. 
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4.2 The Vlasov-Poisson system, the wakefield model and singular solutions 

Besides integrability of the Vlasov-Benney equation, there are other important apphcations of 
the Vlasov equation that have in common the presence of a quadratic term in Aq within the 
Hamiltonian: 

H = ^ jA2{q)dq+]^ j j A^{q) G{q,q') A^{q') dqdq' . 

For example, when G = {dg) ^ , this Hamiltonian leads to the Vlasov-Poisson system, which is 
of fundamental importance in many areas of plasma physics. More generally, this Hamiltonian 
is widely used for beam dynamics in particle accelerators: in this case G is related to the 
electromagnetic interaction of a beam with the vacuum chamber. The wake field is originated 
by the image charges induced on the walls by the passage of a moving particle: while the beam 
passes, the charges in the walls are attracted towards the inner surfaces and generate a field 
that acts back on the beam. This affects the dynamics of the beam, thereby causing several 
problems such as beam energy spread and instabilities. In the literature, the wake function 
W is introduced so that [32] 

G{q,q')= r W{x,q')dx 



Wake functions usually depend only on the properties of the accelerator chamber. 

An interesting wakefield model has been presented in [31] where G is chosen to be the 
Green's function of the Helmholtz operator (l — a^dg): this generates a Vlasov-Helmholtz (VH) 
equation [6] that is particularly interesting for future work. Connections of this equation with 
the well known integrable KdV equation have been proposed. However we believe that this is 
not a natural step since integrability appears already with no further approximations in the 
Vlasov-Benney system that governs the collective motion of the beam. In particular we would 
like to understand the VH equation as a special deformation of the integrable Vlasov-Benney 
case that allows the existence of singular solutions. Indeed, the presence of the Green's 
function G above is a key ingredient for the existence of the single-particle solution, which is not 
allowed in the VB case. In particular, the single-particle solutions for the Vlasov-Helmholtz 
equation may be of great interest, since these singular solutions arise from a deformation of an 
integrable system. The Vlasov-Helmholtz solutions may differ considerably from the well known 
particle behavior of the Vlasov-Poisson case. However, the limit as the deformation parameter 
(the length-scale a) in the Helmholtz Greens function passes to zero (a — > 0) in the wake-field 
equations, one recovers the integrable Vlasov-Benney case. 

4.3 The EPDiff equation and singular solutions 

Another interesting moment equation is given by the integrable EPDiff equation [5j. In this 
case, the Hamiltonian is purely quadratic in the first moments: 

H = \ jjAi{q)G{q,q')AM)dqdq' 

and the EPDiff equation |16| 

^ + ^ 1 G{q,q')A^{q',t)dq' + 2^1 ^ / G{q,q')A^{(^ ,t)dq' = 

comes from the closure of the KMLP bracket given by cotangent lifts. (Without this restriction 
we would obtain again the equations (13. Ih with (3 = G * Ai.) Thus this EPDiff equation is a 



Gibbons, Holm & Tronci Vlasov moments, integrable systems and singular solutions 11 



geodesic equation on the group of diffeomorphisms. The Camassa-Holm equation is a particular 
case in which G is the Green's function of the Helmholtz operator 1 — a^9^. Both the CH 
and the EPDiff equations are completely integrable and have a large number of applications in 
fluid dynamics (shallow water theory, averaged fluid models, etc.) and imaging techniques [19] 
(medical imaging, contour dynamics, etc.). 

Besides its complete integrability the EPDiff equation has the important feature of allowing 
singular delta-function solutions. The connection between the CH and EPDiff equations and the 
moment dynamics lies in the fact that singular solutions appear in both contexts. The existence 
of this kind of solution for EPDiff leads to investigate its origin in the context of Vlasov moments. 
More particularly we wonder whether there is a natural extension of the EPDiff equation to all 
the moments. This would again be a geodesic (hierarchy of) equation, which would perhaps 
explain how the singular solutions for EPDiff arise in this larger context. 



5 Geodesic motion and singular solutions 
5.1 Quadratic Hamiltonians 

The previous examples show how quadratic terms in the Hamiltonian produce interesting be- 
haviour in various contexts. This suggests that a deeper analysis of the role of quadratic terms 
may be worthwhile particularly in connections between Vlasov p— moment dynamics and the 
EPDiff equation, with its singular solutions. Purely quadratic Hamiltonians are considered in 
|15j . leading to the problem of geodesic motion on the space of p— moments. 

we are interested in the problem of geodesic motion on the space of p— moments. In this 
problem the Hamiltonian is the norm on the p— moment given by the following metric and inner 
product. 



l\\Af = II An{q)Gns{q,q')Mq')dqdq' (5.1) 



h 

2 

n,s=0 

The metric Gns{q, q ') in (|5.ip is chosen to be positive definite, so it defines a norm for {A} G 0*. 
The corresponding geodesic equation with respect to this norm is found as in the previous section 
to be, 

= {Am , h} = [nPn-g-Am+n-i + {m + n) — /?„ j (5.2) 

n=0 

with dual variables /Jn G defined by 

/5- = 71-=E / Gnsiq,q')Asiq')dq' = Y,Gns*As. (5.3) 

Thus, evolution under (j5.2p may be rewritten as coadjoint motion on 0* 

dAm ^ 

—^ = {Am,h}=:-22^'ipnAm+n-i (5.4) 

n=0 

This system comprises an infinite system of nonlinear, nonlocal, coupled evolutionary equations 
for the p— moments. In this system, evolution of the m*'^ moment is governed by the potentially 
infinite sum of contributions of the velocities /3„ associated with n*^ moment sweeping the 
(m + n — moment by a type of coadjoint action. Moreover, by equation ()5.3p . each of the 
Pn potentially depends nonlocally on all of the moments. 
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5.2 A geodesic Vlasov equation 

Importantly, geodesic motion for the p— moments is equivalent to geodesic motion for the Euler- 
Poincare equations on the symplectomorphisms (EPSymp) given by the following Hamiltonian 

H[f] = ^ll f{q,p)g{q,p,q',p')f{q\p')dqdpdq'dp' (5.5) 
The equivalence with EPSymp emerges when the function Q is written as 

g{q,q',P,p')=^p''Gnm {q,q')p'"'. 

n,m 

and the corresponding Vlasov equation reads as 

f ^0 

where {•, •} denotes the canonical Poisson bracket. 

Thus, whenever the metric Q for EPSymp has a Taylor series, its solutions may be expressed 
in terms of the geodesic motion for the p— moments. More particularly the geodesic Vlasov 
equation presented here is nonlocal in both position and momentum and is equivalent to the 
vorticity equation in two-dimensions and for a particular choice of the metric. However this 
equation extends to more dimensions and to any kind of geodesic motion, no matters how the 
metric is expressed explicitly. 

5.3 Singular geodesic solutions 

We have now clarified the geometric meaning of the moment equations and we can therefore 
characterize singular solutions, since the geodesic Vlasov equation (EPSymp) essentially de- 
scribes advection in phase space. Indeed, the geodesic Vlasov equation possesses the single 
particle solution 

f{q,p,t) =Y^d{q-Qj{t))d{p-Pj{t)) 
j 

which is a well known singular solution that is admitted whenever the phase-space density is 
advected along a smooth Hamiltonian vector field. This happens, for example, in the Vlasov- 
Poisson system and in the general wakefield model. On the other hand, these singular solutions 
do not appear in the Vlasov-Benney equation. 

Equation ()5.2p admits singular solutions of the form 

N 

Mq,t) = ^P^{t)6{q-Qj{t)) (5.6) 
i=i 

In order to show this is a solution in one dimension, one checks that these singular solutions 
satisfy a system of partial differential equations in Hamiltonian form, whose Hamiltonian couples 
all the moments 

^ oo Af 
n,s=0 j,k=l 
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Explicitly, one takes the pairing of the coadjoint equation 

n,s 

with a sequence of smooth functions {(pm (q)} £ind finally obtains the equations for Qj and Pj 
in canonical form, 

dQj _ dHN dPj _ SHn 

IT ~ dPj ' ir~~ dQj ■ 

These singular solutions of EPSymp are also solutions of the Euler-Poincare equations on 
the diffeomorphisms (EPDiff). In the latter case, the single-particle solutions reduce to the 
pulson solutions for EPDiff [5] . Thus, the singular pulson solutions of the EPDiff equation arise 
naturally from the single-particle dynamics on phase-space. A similar result also holds in higher 
dimensions [15j . 

Further remarks on singular solutions. Another kind of singular solution for the moments 
may be obtained by considering the cold-plasma solution of the Vlasov equation 

3 

Indeed exchanging the variables q ^ pm. the single particle PDF leads to the following expression 

f{q,p,t) = ^i^j{p,t)S{q - Xj{p,t)) 

j 

which is always a solution of the Vlasov equation because of the symmetry in q and p. This 
leads to the following singular solutions for the moments: 

j 

At this point, if one considers a Hamiltonian depending only on Ai (i.e. one considers the action 
of cotangent lifts of diffeomorphism), then it is possible to drop the p-dependence in the A's and 
thereby recover to the singular solutions previously found for eq. (13. Ih . 

5.4 Examples of simplifying truncations and specializations. 

The problem presented by the coadjoint motion equation ()5.4p for geodesic evolution of p— moments 
under EPDiff needs further simplification. One simplification would be to modify the (doubly) 
infinite set of equations in (15. 4|) by truncating the Poisson bracket to a finite set. These moment 
dynamics may be truncated at any stage by modifying the Lie-algebra in the KMLP bracket to 
vanish for weights m + n — 1 greater than a chosen cut-off value. 

For example, if we truncate the sums to m,n = 0, 1, 2 only, then equation ()5.4p produces the 
coupled system of partial differential equations, 

dg{Ao(3^)-2A^dgp2-2/32dgAi 
AodgPo - 2Aidq(3i - (3idqAi - 3A2dgP2 - 2p2dqA2 
2Aidgf3o - SA2d,pi - /3idgA2 



dAo 

dt 

dAi 

"df 
dM 
dt 
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We specialize to the case that each velocity depends only on its corresponding moment, so 
that /3s = G * j4s, s = 0, 1, 2. If we further specialize by setting Aq and A2 initially to zero, then 
these three equations reduce to the single equation 

= -Pid.Ai- 2Aid,Pi. 

Finally, if we assume that G in the convolution Pi = G * Ai is the Green's function for the 
operator relation 

Ai = {l- a'dl)[5i 

for a constant lengthscale a, then the evolution equation for Ai reduces to the integrable 
Camassa-Holm (CH) equation [5] in the absence of linear dispersion. This is the one-dimensional 
EPDiff equation, which has singular (peakon) solutions. Thus, after these various specializations 
of the EPDiff p— moment equations, one finds the integrable CH peakon equation as a further 
specialization of the coadjoint moment dynamics of equation (j5.4p . 

That such a drastic restriction of the p— moment system still leads to such an interesting 
special case bodes well for future investigations of the EPSymp p— moment equations. Further 
specializations and truncations of these equations will be explored elsewhere. Before closing, 
we mention one or two other open questions about the solution behavior of the p— moments of 
EPSymp. 

6 Open questions for future work 

Emergence of singular solutions. Several open questions remain for future work. The first 
of these is whether the singular solutions found here will emerge spontaneously in EPSymp dy- 
namics from a smooth initial Vlasov PDF. This spontaneous emergence of the singular solutions 
does occur for EPDiff. In fact, integrability of EPDiff in one dimension by the inverse scatter- 
ing transform shows that only the singular solutions (peakons) are allowed to emerge from any 
confined initial distribution in that case [5j (this also happens in higher dimensions as it is shown 
by numerical simulations). In contrast, the point vortex solutions of Euler's fluid equations 
(which are isomorphic to the cold plasma singular solutions of the Vlasov-Poisson equation) 
while comprising an invariant manifold of singular solutions, do not spontaneously emerge from 
smooth initial conditions. Nonetheless, something quite analogous to the singular solutions is 
seen experimentally for cold plasma in a Malmberg-Penning trap [llj. Therefore, one may ask 
which outcome will prevail for the singular solutions of EPSymp. Will they emerge from a 
confined smooth initial distribution, or will they only exist as an invariant manifold for special 
initial conditions? One might argue that in two dimensions, the EPSymp equation encompasses 
the equation of vorticity and thus spontaneous emergence of point vortices should not occur. 
However it is possible that the choice of the metric plays an important role in this matter. Of 
course, the interactions of these singular solutions for various metrics and the properties of their 
collective dynamics is a question for future work. 

Possible connections with the Bloch-Iserles system The EPSymp equation is surpris- 
ingly similar in construction to another important integrable geodesic equation on the linear 
Hamiltonian vector fields (Hamiltonian matrices), which has recently been proposed [3J. This 
is a finite dimensional equation whose dynamical variables are symmetric matrices. Now it has 
been shown that this system may be written as the geodesic equation on the group of the linear 
canonical transformations Sp(M; 2n) [4]. This association to canonical transformations raises the 
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question wether it is possible to establish connections with the geodesic Vlasov equation that 
was introduced here 
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